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Abstract 

Background: Low 25-hydroxyvitamin D levels are a global health problem with northern countries such as Canada 
at particular risk. A number of sociodemographic factors have been reported to be associated with low vitamin D 
levels but prior studies have been limited by the ability of the researchers to gather this data directly from clinical 
trial participants. The purpose of this study was to use a novel methodology of inferring sociodemographic 
variables to evaluate the correlates of vitamin D levels in individuals dwelling in the City of Calgary, Alberta, Canada. 

Methods: We utilized data on vitamin D test results from Calgary Laboratory Services between January 1 2010 and 
August 31 201 1. In addition to vitamin D level, we recorded age, sex, and vitamin D testing month as individual-level 
variables. We inferred sociodemographic variables by associating results with census dissemination areas and using 
Census Canada data to determine immigration status, education, median household income and first nations status 
as clustered variables. Associations between vitamin D status and the individual- and dissemination area-specific 
variables were examined using the population-averaged regression model by a generalized estimating equations 
approach to account for the clustering in the data. 

Results: 158,327 individuals were included. Age, sex, month of vitamin D testing (at an individual level), and 
education, immigrant status, first nations status and income (at an aggregate level) were all statistically significant 
predictors of vitamin D status. 

Conclusions: Vitamin D status was associated with a number of sociodemographic variables. Knowledge of these 
variables may improve targeted education and public health initiatives. 
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Background 

Vitamin D has received considerable interest from the 
medical community and the public at large because of its 
importance in the maintenance of health [1-6] combined 
with the finding of widespread global deficiency [7-12]. 
Among the major micronutrients, vitamin D is unique in 
that the primary source is solar ultraviolet B radiation 
[13] (which converts 7-dehydocholesterol to vitamin D 3 
in the skin), with oral supplementation representing an 
important secondary source [14]. These characteristics 
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make vitamin D deficiency an important target for public 
health interventions. 

High latitude countries such as Canada may be espe- 
cially vulnerable to vitamin D deficiency because of 
lower ultraviolet radiation levels [13-23]. Indeed findings 
of a Canadian cross sectional survey suggested that as 
many as three million Canadians have inadequate vita- 
min D levels and 1.1 million Canadians are vitamin D 
deficient (defined as a serum vitamin D level <27.5 nmol/L; 
to convert 25-hydroxyvitamin D levels from SI to con- 
ventional units divide by 2.496) [17]. Other specific 
sociodemographic factors have been reported to be asso- 
ciated with low vitamin D levels including advanced age 
and/or residence in a nursing home [24], cultural factors 
such as skin covering [25,26], aboriginal ancestry [27], low 
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dietary intake of milk [17], darker skin pigmentation and/ 
or non-white ethnicity [14,17,18,21,28], obesity [18,29], and 
lower education [30] . 

However there is scope for additional large scale 
examination of the sociodemographic correlates of vita- 
min D status, particularly with regards to aboriginal an- 
cestry, education and the largely unstudied variable of 
household income, as these variables have received only 
limited attention in prior studies [28,30]. In this study 
we combined a secondary analysis of laboratory test re- 
sults with aggregate census Canada data to determine 
sociodemographic factors independently associated with 
serum vitamin D levels in a large sample of individuals 
in Calgary, a northern Canadian city of Alberta province. 
Although prior researchers have utilized geospatial map- 
ping to infer health-related variables [31-34], this ap- 
proach is novel in that it is being applied for the first 
time to vitamin D data. We hypothesized that spatial 
variance in sociodemographic factors within the city of 
Calgary will be associated with spatial variance in mean 
25-hydroxyvitamin D levels, and that these variances 
can be used to infer sociodemographic associations with 
25-hydroxyvitamin D level. 

The recent Institute of Medicine (IOM) report on 
Dietary Reference Intakes for vitamin D for Canada and 
the USA, suggests that for the skeletal benefits of vita- 
min D, a 25-hydroxyvitamin D level level of 50 nmol/L 
reflects adequate vitamin D intake for 97.5% of the 
population [35]. In our analysis, we therefore defined 
25-hydroxyvitamin D sufficiency as a serum level of 
greater then or equal to 50 nmol/L. 

Methods 

Ethics statement 

The study protocol was approved by the University of 
Calgary Conjoint Health Review Ethics Board (Ethics 
ID 23919). 

Study population and data sources 

We undertook this observational study combining la- 
boratory data with clustered (census dissemination area 
level) variables obtained from the 2006 Canadian Census 
of Population. The study population consisted of adults 
25 years of age and older who underwent vitamin D testing 
at Calgary Laboratory Services (CLS) between January 01, 
2010 and August 31, 2011. CLS is the sole provider of 
laboratory testing to Calgary, Alberta and the surround- 
ing areas (approximate population 1.4 million). All 25- 
hydroxyvitamin D tests were performed as part of routine 
patient care and were analyzed in a single laboratory using 
the LIASON 25-hydroxyvitamin D Total assay (Diasorin, 
Ltd.) on Roche modular analyzers. The lower limit of de- 
tection of this assay was 10 nmol/L. For the purpose of 
analysis, values <10 nmol/L were recorded as 10 nmol/L. 



Quality control is performed daily on these analyzers. 
Quality assurance is performed through subscription to the 
Vitamin D External Quality Assessment Scheme (DEQAS). 
For the year 2012 our lab's average bias from our method 
mean (DiaSorin Liaison Total) for all samples was less than 
5%, with no individual result being more than 10% from 
the method mean. The inter-assay coefficient of variation 
was 4.2% at 42 nmol/1 and 2.9% at 130 nmol/L. Intra- 
Assay coefficient of variations from our daily quality con- 
trols was 8.5% at 50 nmol/L, 8.3% at 85 nmol/L and 9.6% 
at 265 nmol/L. 

To avoid pseudoreplication, each patient was included 
only once in the analysis. If more than one record 
existed for the same individual, one record was chosen 
at random. An alternative approach would have been to 
choose the first available test result for each individual; 
this may have excluded some measurements taken to de- 
termine response to supplementation. However, because 
of the limited time frame covered by our data and our 
inability to determine if patients were taking supple- 
ments, we opted to use one test result taken at random 
for each patient. For each vitamin D test record the fol- 
lowing information was extracted from the laboratory 
information system: 25-hydroxyvitamin D level, age, sex, 
month of testing, and health care number. Health care 
number was then used as a linking variable to subject 
postal codes from an Alberta Health Services database. 
The postal codes were converted to their corresponding 
census dissemination area and geographic coordinates 
using a Statistics Canada Postal Code Conversion File. 
Census Dissemination Areas (CDA) are the lowest geo- 
graphic level in the Canadian Census and consist of 
geographical groupings of 300-700 people. We consid- 
ered only individuals residing within the City of Calgary 
and so census dissemination areas outside of the city 
limits were removed from the dataset. Finally, 25- 
hydroxyvitamin D levels greater than 300 nmol/L were 
deemed as extreme values and removed. 

The resulting dataset was plotted onto a map of the 
city of Calgary using ArcGIS v. 9.3 software at the 
Spatial and Numeric Services Department of the Uni- 
versity of Calgary Library. Following mapping of the 
individual records, sociodemographic variables for each 
dissemination area within the city of Calgary were 
imported from Statistics Canada data obtained from the 
2006 Canadian Census [36]. These data were then joined 
as clustered variables at the level of census dissemination 
area from the 2006 Canada Census. The following dis- 
semination area variables were added to the dataset: 
percent of individuals of aboriginal descent, median 
household income, percent of individuals born outside 
of Canada, and percent of individuals over age 25 
with at least some post-secondary education. This com- 
bined data set was then used to examine associations 
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between 25-hydroxyvitamin D level and the various 
sociodemographic variables. 

Statistical analysis 

Data were presented as mean and standard deviation 
(SD), or median and interquartile range (IQR), when ap- 
propriate, for continuous variables; and, as frequencies 
and proportions for categorical variables. Given the cor- 
related (clustered) nature of the data (i.e. individuals 
were clustered within census dissemination areas where 
they lived in as the data was collected), we used the SAS 
macro % RCS_Reg as detailed in [37] to fit the marginal 
(population-averaged) generalized estimating equations 
(GEE) regression model to account for the clustering in 
the data, while using restricted cubic spline (RCS) func- 
tions (i) to visually and statistically check the assumption 
of linearity of the association between each of the con- 
tinuous predictors and the outcome, and, (ii) to graphic- 
ally characterize and quantify the association when the 
latter assumption is not valid [37]. The number of knots 
for the splines of each continuous predictor was chosen 
according to the information criteria (QIC and QICu, 
the smaller the better). The final GEE model included 
each of the continuous predictors with their restricted 
cubic splines, as well as gender and month of vitamin D 
testing as categorical variables. All reported P values 
were two-sided, and considered as significant if < 0.05. 

Results and discussion 

Data were abstracted from a total of 1986 census dis- 
semination areas within the City of Calgary, with a me- 
dian number of 111 (IQR: 68-243) individuals being 
surveyed in each area. Outlier 25-hydroxyvitamin D 
levels greater than 300 were excluded from the analysis 
(119 cases) resulting in individual level data for a total of 
158,327 patients. The majority (64.9%) were female, and 
most (73.9%) had a vitamin D level greater than or equal 
to 50 nmol/L. Median 25-hydroxyvitamin D levels as 
well as the proportion of individuals falling below the In- 
stitute of Medicine categories of vitamin adequacy [35] 
are given in Table 1. Associations of 25-hydroxyvitamin 

Table 1 Median 25 hydroxyvitamin D level for the Calgary 
population studied, and the proportion of the population 
below levels representing the Institute Of Medicine's [31] 
categories of vitamin adequacy: deficiency (<30 nmol/L), 
estimated average requirement (40 nmol/L) and 



recommended dietary allowance (50 nmol/L) 





Males 


Females 


Median 25 hydroxyvitamin D level 


65 


71 


Proportion of individuals <30 nmol/L 


0.10 


0.07 


Proportion of individuals <40 nmol/L 


0.19 


0.15 


Proportion of individuals <50 nmol/L 


0.30 


0.24 



area-level variables are summarized in Table 2. The 
mean vitamin D level per dissemination area is shown in 
Figure 1. As this map shows, higher average levels of 
vitamin D tended to be found in the inner city neighbor- 
hoods while the lowest levels tended to cluster in the far 
north and northeast. Chloropleths of median household 
income, percent of the population with at least some 
post-secondary education, percent of individuals of 
aboriginal decent and percent of individuals born out- 
side of Canada (Figure 2) revealed some broad socio- 
demographic trends within the city. Notably, higher 
incomes and higher education levels tended to cluster 
in inner city neighborhoods and in the northwest (near 
the University of Calgary), while immigrants tended to 
be found in the north and northeast. The distribution 
of aboriginal people tended to be highest in the east 
and in certain inner city neighborhoods. 

Table 3 shows the results from the population-averaged 
GEE model investigating the association between 25- 
hydroxy vitamin D status and each of the predictors. The 
test for the overall association between 25-hydroxyvitamin 
D level and each of the continuous predictors was signifi- 
cant (all P-values < 0.0001). The test for non-linear associ- 
ation was also significant, which means that the association 
between 25-hydroxyvitamin D level and each of the con- 
tinuous predictors was significantly not linear (p-value for 
linearity test was <0.0001, 0.0295, <0.0001, 0.009 and 
0.0002 for age, median household income, percent of indi- 
viduals born outside of Canada, percent of individuals of 
aboriginal decent and percent of population with at least 
some post-secondary education, respectively). Figure 3 
graphically displays the nonlinear relationships between 
25-hydroxyvitamin D level and each of the continuous 
predictors, with each graph showing the adjusted differ- 
ence in mean 25-hydroxyvitamin D level (as indicated 
on the y-axis) between individuals with any value of the 
corresponding continuous predictor (as indicated on the 
x-axis) and individuals with the specified reference value 
(as described in Figure 3), when the other continuous vari- 
ables plus gender and month of 25-hydroxyvitamin D test- 
ing were adjusted for. Interestingly, as shown in Figure 3, 
when other covariates were controlled in the generalized 
estimating equation model, median household income and 
percent of individuals of aboriginal descent showed only 
very modest associations with 25-hydroxyvitamin D level, 
although the associations were statistically significant. In 
contrast, education showed a clear relationship in that the 
mean levels of 25-hydroxyvitamin D were higher in census 
dissemination areas that had a greater proportion of indi- 
viduals with at least some post-secondary education. Simi- 
larly, census dissemination area immigrant status was 
related to lower census dissemination area levels of 25- 
hydroxyvitamin D. The function for age was more complex 
with a nadir at age around 35 years then a steady increase 



Naugler et al. BMC Public Health 2013, 13:316 
http://www.biomedcentral.eom/1 471-2458/1 3/3 1 6 



Page 4 of 10 



Table 2 Data summary of the individual- and census dissemination area-level variables included in the analysis, 
stratified by the presence of vitamin D insufficiency 



Variable 


Vitamin D insufficient (< 50 nmol/L) 


Vitamin D sufficient (> 50 nmol/L) 


Total 




(N=41,401) 




(N=1 16,926) 


(N=1 58327) 


Age (years) 


Mean (SD) 


47.7 (14.4) 




53.6 (15.4) 


52.1 (15.4) 


Gender, n (%) 


Female 


24262 (58.6) 




78417 (67.1) 


102679 (64.9) 


Month of Vitamin D testing, n (%) 


January 


4084 (9.9) 




10615 (9.1) 


14699 (9.3) 


Feburary 


4210 (10.2) 




10996 (9.4) 


15206 (9.6) 


March 


6025 (14.6) 




12657 (10.8) 


18682 (11.8) 


April 


4371 (10.6) 




11153 (9.5) 


15524 (9.8) 


May 


3999 (9.7) 




12245 (10.5) 


16244 (10.3) 


June 


4345 (10.5) 




12409 (10.6) 


16754 (10.6) 


July 


2939 (7.1) 




10736 (9.2) 


13675 (8.6) 


August 


2863 (6.9) 




10912 (9.3) 


13775 (8.7) 


September 


2015 (4.9) 




6723 (5.7) 


8738 (5.5) 


October 


2265 (5.5) 




6745 (5.8) 


9010 (5.7) 


November 


2277 (5.5) 




6618 (5.7) 


8895 (5.6) 


December 


2008 (4.9) 




5117 (4.4) 


7125 (4.5) 


Median household income of the dissemination area ($CDN) 


Median (IQR 1 ) 


70690 (52623-87577) 




73831 (55649-96393) 


73264 (54982-94482) 


Percent of individuals in the dissemination area not born in Canada 


Median (IQR) 


25.6 (17.3-37.1) 




21.0 (14.1-29.9) 


22.3 (14.8-31.9) 


Percent of aboriginals in the dissemination area 


Median (IQR) 


1 .7 (0.0 - 3.4) 




1 .6 (0.0 - 3.0) 


1.6 (0.0-3.1) 


Percent of individuals >= 25 years old in the dissemination area with an university degree 


Median (IQR) 


49.0 (35.9-62.6) 




56.9 (43.9-67.9) 


55.7 (41.1 - 66.7) 



1 IQR=interquartile range. 



in vitamin D level before leveling off around age 65 years 
(Figure 3). 

In addition, gender had also a statistically significant 
association with 25-hydroxyvitamin D level, with fe- 
males exhibiting on average a 25-hydroxyvitamin D 
level 7.5 nmol/L (95% CI: 7.1 - 7.9 nmol/L, p < 0.0001) 
higher than males. We observed that mean vitamin D 
level differed according to month (Table 1), although 
the investigation of these differences was not one of 
purposes of this paper. 

Overall, 26% of individuals in our study were 25- 
hydroxyvitamin D deficient or insufficient (serum 25- 
hydroxyvitamin D < 50 nmol/L), compared with 20% of 
individuals in the Canadian multicentre osteoporosis 
study (CAMOS) [18]. The difference may be due to the 
differing ages of the subjects in these studies: whereas 
the (CAMOS) only considered adults over age 35, we in- 
cluded individuals age 25 and older. However, it should 
also be noted that the CAMOS subjects were recruited 
by random selection from the population living within 



50 km of the centre of the Canadian cities studied, while 
the present study is not a random selection of the 
Calgary population. It is also possible that the difference 
may be due to different proportions of non-white sub- 
jects in each study. The 2007-2009 Canadian Health 
Measures Survey [17] reported a relationship between 
age and 25-hydroxyvitamin D level with the 20-39 
age groups for both women and men showing the 
lowest vitamin D levels. This is also consistent with 
our findings. 

Our observation of a nadir in 25- OH vitamin D 
levels in young adults may seem at odds with the 
widely held perception that vitamin D deficiency is 
principally a problem of older individuals [24,38-41]. It 
is possible that the low levels of 25-hydroxyvitamin D 
in the elderly may not be related to age per se but to 
other correlated variables such as education. Alterna- 
tively, it is possible that older individuals or those with 
greater post-secondary education are achieving higher 
25-hydroxyvitamin D levels through supplementation. 
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The inferred relationship we found between immi- 
grant status and lower 25-hydroxyvitamin D levels is 
likely due to immigrant status being a proxy measure 
for increased skin pigmentation, a well-known risk factor 
for 25-hydroxyvitamin D deficiency [17,19,21,42-44]. The 
previously reported lower 25-hydroxyvitamin D level in 
Canadian aboriginals [27] is only weakly supported by 
our data. Rather, our findings suggest that the reported 
association may be due to other correlated factors such 
as age and education level. We confirmed this by com- 
paring the crude association between the percentage of 



aboriginal people and vitamin D levels with the adjusted 
association between percentage of aboriginal people and 
vitamin D levels. The resulting crude association was 
stronger (i.e., showed a steeper slope between percentage 
of aboriginal people and vitamin D levels) than the ad- 
justed association, suggesting that the association observed 
between percentage of aboriginal people and vitamin D 
could be indeed explained by age and education level (ana- 
lysis not shown). 

Interestingly, income demonstrated a very minor 
association with 25-hydroxyvitamin D level. This was 
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Table 3 Estimates from the population-averaged GEE model investigating associations with vitamin D level 



Parameter 


Estimate 


Standard error 


P-value 


Overall p-value 


Intercept 


64.5883 


2.4201 


<.0001 


<.0001 


Age (years) 








<.0001 


age_RCS_lin 


-0.2738 


0.0422 


<.0001 




age_RCS_S1 


0.0015 


0.0001 


<.0001 




age_RCS_S2 


-0.0030 


0.0004 


<.0001 




age_RCS_S3 


0.0008 


0.0004 


0.0823 




Median household income ($CDN 10000) 








<.0001 


income_RCS_lin 


-0.4078 


0.2350 


0.0826 




income_RCS_S1 


0.0543 


0.0263 


0.0392 




income_RCS_S2 


-0.1270 


0.0793 


0.1092 




income_RCS_S3 


0.0824 


0.0855 


0.3352 




Percentage of immigrants 








<.0001 


percentimmi RCSJin 


-7.6190 


9.4844 


0.4218 




percentimmi RCS_S1 


-172.823 


567.4589 


0.7607 




percentimmi RCS_S2 


-1 122.79 


1547.330 


0.4681 




percentimmi RCS_S3 


3152.879 


151 1.398 


0.0370 




Percentage of aboriginal people 








<.0001 


percentabor RCSJin 


-68.5534 


15.3580 


<.0001 




percentabor RCS_S1 


13868.36 


531 1.648 


0.0090 




Percentage of subjects with some post-secondary education 








<.0001 


percentuniv_RCS_lin 


13.6596 


4.2690 


0.0014 




percentuniv_RCS_S1 


237.8126 


69.0070 


0.0006 




percentuniv_RCS_S2 


-842.990 


263.7487 


0.0014 




percentuniv_RCS_S3 


1 090.294 


509.4453 


0.0323 




Gender 








<.0001 


Female 


7.4674 


0.1930 


<.0001 




Male 


0.0000 


0.0000 






Month of vitamin D testing 








<.0001 


December 


-1.3870 


0.5024 


0.0058 




November 


0.1357 


0.4583 


0.7672 




October 


-1.3138 


0.4593 


0.0042 




September 


-0.1319 


0.4590 


0.7738 




August 


4.7383 


0.4252 


<.0001 




July 


3.7308 


0.4142 


<.0001 




June 


0.5034 


0.3959 


0.2035 




May 


3.5073 


0.4330 


<.0001 




April 


0.6403 


0.4129 


0.1209 




March 


-3.6969 


0.4045 


<.0001 




February 


0.4483 


0.4108 


0.2751 




January 


0.0000 


0.0000 







Note: Age, median household income, percentage of immigrants, and percentage of subjects with some post-secondary education were coded using restricted 
cubic spline (RCS) functions with five knots located at the 5th, 25th, 50th, 75th, and 95th percentiles of their distribution, and percentage of aboriginal people was 
coded using the 5th, 50th and 95th percentiles of its distribution. The variables listed underneath each of the continuous variables are their corresponding spline 
variables (e.g., age_RCS_lin, age_RCS_S1, age_RCS_S2, and age_RCS_S3 are the 4 spline variables of age). 
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Figure 3 Adjusted associations of vitamin D level (nmol/L). (A) age (reference value = 37 years); (B) median household income (reference 
value = $CDN 73,264, the median); (C) percentage of individuals who were not born in Canada (reference value = 22%, the median); (D) percentage of 
Aboriginal people (reference value = 2%, the median); (E) percentage of individuals >25 years old who had at least some post-secondary education 
(reference value = 56%, the median). Age, median household income, percentage of immigrants, and percentage of subjects with some post-secondary 
education were coded using restricted cubic spline (RCS) functions with five knots located at the 5th, 25th, 50th, 75th, and 95th percentiles of their 
distribution, and percentage of aboriginal people was coded using the 5th, 50th and 95th percentiles of its distribution. Y-axis represents the difference 
in mean vitamin D level between individuals with any value of the continuous predictor (as indicated on the x-axis) and individuals with the reference 
value (as indicated by the vertical dashed line). Shaded area represents the 95% confidence band. Horizontal dashed line indicates the zero difference in 
mean 25 hydroxyvitamin D level. 



somewhat unexpected as income has been shown to be 
an important health predictor in other contexts [45]. 
Finally, the association between 2 5 -hydroxyvitamin D 
level and post-secondary education is consistent with a 
recent Finnish study [30]. It is tempting to speculate 
that this relationship may be related to increased use of 
vitamin D supplements among individuals with higher 
education [46]. 

Our study should be interpreted in light of its limitations. 
Firstly, because we relied on the secondary use of labora- 
tory results combined with Census Canada data, we were 
unable to capture potentially important information such 
as dietary milk consumption [17,18,47], body mass 
index [18,29] and 25 -hydroxyvitamin D supplementa- 
tion [18]. Nor did we have access to a surrogate for 
25-hydroxyvitamin D status such as parathyroid hormone 



level. A second potential limitation concerns the changing 
nature of the demographic make-up of the city of Calgary. 
Although we used 25-hydroxyvitamin D results from 2010 
to 2011, the census Canada data we used was from the 
2006 census. If there was a considerable demographic shift 
within a given census dissemination area in the intervening 
few years, this could have affected the clustered variables 
we used. Third, although dramatic differences in mean 
25-hydroxyvitamin D levels existed among census dis- 
semination areas, we cannot exclude the possibility of a 
confounding effect if the probability of individuals with low 
or high vitamin D being tested also varied with these 
socio-demographic variables. However, our large sam- 
ple size should mitigate against this effect. Finally, this 
study shares the potential weakness of all studies util- 
izing ecological data in that the inferences regarding 
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group level variables may not necessarily reflect indi- 
vidual level variables. 

Conclusions 

In this study we examined the associations between a num- 
ber of sociodemographic variables and 25-hydroxyvitamin 
D level using a combination of secondary clinical data with 
individual level variables and clustered variables derived 
from Census Canada data. We found that the mean levels 
of 25-hydroxyvitamin D for the city of Calgary varied 
widely by census dissemination area and that the predom- 
inant predictors of this variation seemed to be age, educa- 
tion level and immigration status among the variables 
considered in this study. 

Competing interests 

The authors declare no relevant competing interests. 
Authors' contributions 

CN conceived of the study. All authors contributed to the study design. CN 
collected the data. CN and JZ performed he analyses. CN drafted the 
manuscript. All authors contributed to revisions, read and approved the final 
manuscript. 

Acknowledgements 

The authors wish to thank Mr. Peter Peller, Spatial and Numeric Services, 
University of Calgary Library for assistance with the geospatial analysis, and 
Megan-Joy Rockey for help with data management. 

Author details 

department of Pathology and Laboratory Medicine, University of Calgary, 
Calgary, Alberta, Canada, department of Family Medicine, University of 
Calgary, Calgary, Alberta, Canada. 3 Calgary Laboratory Services, Calgary, 
Alberta, Canada. 4 C414, Diagnostic and Scientific Centre, 9, 3535 Research 
Road NW, Calgary AB T2L 2K8, Canada, department of Medicine, University 
of Calgary, Calgary, Alberta, Canada. 

Received: 15 October 2012 Accepted: 4 April 2013 
Published: 8 April 2013 

References 

1 . Holick MF: Vitamin D deficiency. N Engl J Med 2007, 357:266-281 . 

2. Bischoff-Ferrari HA, Willett WC, Wong JB, Stuck AE, Staehelin HB, Orav EJ, 
Thoma A, Kiel DP, Henschkowski J: Prevention of nonvertebral fractures 
with oral vitamin D and dose dependency: a meta-analysis of 
randomized controlled trials. Arch Intern Med 2009, 169:551-561. 

3. Bischoff-Ferrari HA, Dawson-Hughes B, Staehelin HB, Orav JE, Stuck AE, 
Theiler R, Wong JB, Egli A, Kiel DP, Henschkowski J: Fall prevention with 
supplemental and active forms of vitamin D: a meta-analysis of 
randomised controlled trials. Br Med J 2009, 339:b3692. 

4. Lin J, Manson JE, Lee IM, Cook NR, Buring JE, Zhang SM: Intakes of calcium 
and vitamin D and breast cancer risk in women. Arch Intern Med 2007, 
167:1050-1059. 

5. Gorham ED, Garland CF, Garland FC, Grant WB, Mohr SB, Lipkin M, Newmark HL, 
Giovannucci E, Wei M, Holick MF: Optimal vitamin D status for colorectal 
cancer prevention: a quantitative meta analysis. Am J Prev Med 2007, 
32:210-216. 

6. Giovannucci E, Liu Y, Hollis BW, Rimm EB: 25-hydroxyvitamin D and risk of 
myocardial infarction in men: a prospective study. Arch Intern Med 2008, 

168:1174-1180. 

7. Moan J, Porojnicu AC, Dahlback A, Setlow RB: Addressing the health 
benefits and risks, involving vitamin D or skin cancer, of increased sun 
exposure. PNAS 2008, 105:668-673. 

8. Pieper CF, Colon-Emeric C, Caminis J, Betchyk K, Zhang J, Janning C, 
Shostak J, LeBoff MS, Heaney RR, Lyles KW: Distribution and correlates of 
serum 25-hydroxyvitamin D levels in a sample of patients with hip 
fracture. Am J Geriatr Pharmacother 2007, 5:335-340. 



9. Venning G: Recent developments in vitamin D decifiency and muscle 
weakness among elderly people. Br Med J 2005, 330:524-526. 

1 0. Meunier P: Prevention of hip fractures by correcting calcium and vitamin D 
insufficiencies in elderly people. Scand J Rheumatol 1996, 25(Suppl 103):75-78. 

1 1 . Hayes CE: Vitamin D: a natural inhibitor of multiple sclerosis. Proc Nurt 
Soc 2000, 59:531-535. 

12. Bischoff-Ferrari HA, Can U, Staehelin HB, Platz A, Henschkowski J, Michel BA, 
Dawson-Hughes B, Theiler R: Severe vitamin D deficiency in Swiss hip 
fracture patients. Bone 2008, 42:597-602. 

13. Kimlin MG: Geographic location and vitamin D synthesis. Mol Aspects Med 

2008, 29:453-461. 

14. Chen TC, Chimeh F, Lu Z, Mathieu J, Person KS, Zhang A, Kohn N, Martinello S, 
Berkowitz R, Holick MF: Factors that influence the cutaneous synthesis and 
dietary sources of vitamin D. Arch Biochem Biophys 2007, 460:21 3-21 7. 

15. Touma Z, Eder L, Zisman D, Feld J, Chandran V, Rosen CF, Shen H, Cook RJ, 
Gladman DD: Seasonal variation in vitamin D levels in psoriatic arthritis 
patients from different latitudes and its association with clinical 
outcomes. Arthritis Care Res (Hoboken) 201 1, 63:1440-1447. 

16. Rucker D, Allan JA, Fick GH, Hanley DA: Vitamin D insufficiency in a 
population of healthy western Canadians. CMAJ 2002, 166:1517-1524. 

17. Langolois K, Greene-Finestone L, Little J, Hidiroglou N, Whiting S: Vitamin D 
status of Canadians as measured in the 2007 to 2009 Canadian Health 
Measures Survey. Statistics Canada 2010, 21:47-54. 

18. Greene-Finestone LS, Berger C, de Groh M, Hanley DA, Hidiroglou N, Sarafin K, 
Poliquin S, Krieger J, Richards JB, Goltzman D, CaMos Research Group: 
25-Hydroxyvitamin D in Canadian adults: biological, environmental, 
and behavioral correlates. Osteoporos Int 2011, 22:1389-1399. 

19. Vieth R, Cole DE, Hawker GA, Trang HM, Rubin LA: Wintertime vitamin D 
insufficiency is common in young Canadian women and their vitamin D 
intake does not prevent it. Eur J Clin Nutr 2001, 55:1091-1097. 

20. Mark S, Gray-Donald K, Delvin EE, O'Loughlin J, Paradis G, Levy E, Lambert M: 
Low vitamin D status in a representative sample of youth from Quebec, 
Canada. Clin Chem 2008, 54:1283-1289. 

21. Gozdzik A, Barta JL, Wu H, Wagner D, Cole DE, Vieth R, Whiting S, Parra 
EJ: Low wintertime vitamin D levels in a sample of healthy young 
adults of diverse ancestry living in the Toronto area: Associations 
with vitamin D intake and skin pigmentation. BMC Public Health 2008, 
8:336. 

22. Newhook LA, Sloka S, Grant M, Randell E, Kovacs CS, Twells LK: Vitamin D 
insufficiency common in newborns, children and pregnant women living 
in Newfoundland and Labrador, Canada. Matern Child Nutr 2009, 5:186-191. 

23. Waiters B, Godel JC, Basu TK: Perinatal vitamin D and calcium status of 
northern Canadian mothers and their newborn infants. J Am Coll Nutr 
1999, 18:122-126. 

24. van der Wielen RP, Lowik MR, van den Berg H, de Groot LC, Haller J, 
Moreiras O, van Staveren WA: Serum vitamin D concentrations among 
elderly people in Europe. Lancet 1995, 346:207-210. 

25. Dawodu A, Kochiyil J, Altaye N: Pilot study of sunlight exposure and 
vitamin D status in Arab women of childbearing age. East Mediterr Health 
J 2011, 17:570-574. 

26. Budak N, Cicek B, Sahin H, Tutus A: Bone mineral density and serum 
25-hydroxyvitamin D level: is there any difference according to the 
dressing style of the female university students. Int J Food Sci Nutr 
2004, 55:569-575. 

27. Weiler HA, Leslie WD, Krahn J, Steiman PW, Metge CJ: Canadian Aboriginal 
women have a higher rate of vitamin D deficiency than non-Aboriginal 
women status despite similar dietary vitamin D intakes. J Nutr 2007, 

137:461-465. 

28. Ginde AT, Liu MC, Camargo CA Jr: Demographic differences and trends of 
vitamin D insufficiency in the US population, 1988-2004. Arch Intern Med 

2009, 169:626-632. 

29. Saintonge S, Bang H, Gerber LM: Implications of a new definition of vitamin D 
deficiency in a multiracial US adolescent population: the National Health and 
Nutrition Examination Survey III. Pediatrics 2009, 123:797-803. 

30. Jaaskelainen T, Knekt P, Marniemi J, Sares-Jaske L, Mannisto S, 
Heliovaara M, Jarvinen R: Vitamin D status is associated with 
sociodemographic factors, lifestyle and metabolic health. Eur J Nutr 
2012 [Epub ahead of print]. 

31. Hall ES, South AP, Jones DE, Meinzen-Derr J, Huo S, Liu L, Greenberg JM: 
Spatial analysis in support of community health intervention strategies. 
AMIA Annu Symp Proc 2012, 2012:31 1-320. 



Naugler et al. BMC Public Health 2013, 13:316 
http://www.biomedcentral.eom/1 471-2458/1 3/3 1 6 



Page 10 of 10 



32. Musio M, Sauleau EA: Augustin NH Resources allocation in healthcare for 
cancer: a case study using generalised additive mixed models. Geospat 
Health 2012 Nov, 7(1):83-9. 

33. Rushton G, Krishnamurthy R, Krishnamurti D, Lolonis P, Song H: The spatial 
relationship between infant mortality and birth defect rates in a U.S. 
city. Stat Med 1 996, 1 5(1 7-1 8):1 907-1 9. 

34. DeFranco EA, Lian M, Muglia LA, Schootman M: Area-level poverty and 
preterm birth risk: a population-based multilevel analysis. BMC Public 
Health 2008, 8:316. 

35. IOM (Institute of Medicine): Dietary Reference Intakes for Calcium and Vitamin 
D. Washington, DC: The National Academies Press; 201 1. 

36. Statistics Canada: 2006 Canadian Census. Available at: www.statcan.gc.ca. 

37. Desquilbet L, Mariotti F: Dose-response analyses using restricted cubic 
spline functions in public health research. Stat Med 2010, 29:1037-1057. 

38. Liu BA, Gordon M, Labranche JM, Murray TM, Vieth R, Shear NH: Seasonal 
prevalence of vitamin D deficiency in institutionalized older adults. J Am 
GeriatrSoc 1997, 45:598-603. 

39. Peacock M, Liu G, Carey M, McClintock R, Ambrosius W, Hui S, Johnston CC: 
Effect of calcium or 250H vitamin D3 dietary supplementation on bone 
loss at the hip in men and women over the age of 60. J Clin Endocrinol 
Metab 2000, 85:3011-3019. 

40. Holick MF, Matsuoka LY, Wortsman J: Age, vitamin D, and solar ultraviolet. 
Lancet 1989, 2:1104-1105. 

41. von Muhlen DG, Greendale GA, Garland CF, Wan L, Barrett-Connor E: 
Vitamin D, parathyroid horomone levels and bone mineral density in 
community-dwelling older women: the Rancho Bernardo study. 
Osteoporos Int 2005, 1 6:1 721 -1 726. 

42. Mithal A, Wahl DA, Bonjour JP, Burckhardt P, Dawson-Hughes B, Eisman JA, 
El-Hajj Fuleihan G, Josse RG, Lips P, Morales-Torres J, IOF Committee of 
Scientific Advisors (CSA) Nutrition Working Group: Global vitamin D status and 
determinants of hypovitaminosis D. Osteoporosis Int 2009, 20:1807-1820. 

43. Lebrun JB, Moffatt ME, Mundy RJ, Sangster RK, Postl BD, Dooley JP, Dilling LA, 
Godel JC, Haworth JC: Vitamin D deficiency in a Manitoba community. Can J 
Public Health 1 993, 84:394-396. 

44. Ward LM, Gaboury I, Ladhani M, Zlotkin S: Vitamin D-deficiency rickets 
among children in Canada. CMAJ 2007, 177:161-166. 

45. Lynch JW, Smith GD, Kaplan GA, House JS: Income inequality and mortality: 
importance to health of individual income, psychosocial environment, or 
material conditions. Br Med J 2000, 320:1 200-1 204. 

46. Rock CL: Multivitamin-multimineral supplements: who uses them? Am J 
Clin Nutr 2007, 85:277s-279s. 

47. Vatanparast H, Calvo MS, Green TJ, Whiting SJ: Despite mandatory 
fortification of staple foods, vitamin D intakes of Canadian children and 
adults are inadequate. J Steroid Biochem Mol Biol 2010, 121:301-303. 



doi:1 0.1 1 86/1 471 -2458-1 3-3 1 6 

Cite this article as: Naugler et al.: Association of vitamin D status with 
socio-demographic factors in Calgary, Alberta: an ecological study using 
Census Canada data. BMC Public Health 2013 13:316. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at f~\ RiftMM i rpntral 

www.biomedcentral.com/submit momea central 



